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ABSTRACT 

A generalized hypothesis testing approach is applied to the problem of tracking 
several objects where several different associations of data with objects are pos- 
sible. Such problems occur, for instance, when attempting to distinctly track 
several aircraft maneuvering near each other or when tracking ships at sea. Con- 
ceptually the problem is solved by first associating data with objects in a statis- 
tically reasonable fashion and then tracking with a bank of Kalman filters. 

The objects are assumed to have motion characterized by a fixed but unknown 
deterministic portion plus a random process portion modeled by a shaping filter. 
For example, the object mi^t be assumed to have a mean strai^t line path about 
which it maneuvers in a random manner. Several hypothesized associations of 
data with objects are possible because of ambiguity as to which object the data 
comes from, false alarm/detection errors, and possible uncertainty in the number 
of objects being tracked. 

The statistical likelihood function is computed for each possible hypothesized 
association of data with objects. Then the generalized likelihood is computed 
by maximizing the likelihood over parameters that define the deterministic 
motion of the object. This forms the basis of the generalized hypothesis test- 
ing approach. 

The computational burden is dominated by combinatoric considerations. Proce- 
dures are addressed that relieve the computational burden. 
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1 . Introduction 


This paper deals with the problem of tracking multiple objects 
where it is unclear which measurements can be associated with which 
objects. Such problems arise, for instance, when tracking aircraft 
maneuvering near each other or when tracking ships at sea. The pro- 
blem is approached by first performing data-to-object association in a 
statistically reasonable fashion and then tracking with a bank of 
Kalman filters. 

A recent survey paper Reference 1 describes the various Bayesian 
and likelihood approaches that have been developed to handle multiple 
object tracking problems. This paper follows the likelihood line of 
development pioneered by Sittler in Reference 2 and carried forward by 
various researchers, notably by Morefield in Reference 3. The con- 
tribution of this paper is the use of a generalized likelihood approach 
to deal with unknown system parameters. The conceptual approach is 
applicable to a wide variety of parameters including process noise 
covariances and time constants in shaping filter models for object 
motion. The specific theory developed here applies to unknown 
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parameters that enter measurements linearly. Examples of such para- 
meters are position and velocity coordinates at track initiation or 
parameters describing a fixed but unknown mean path for object motion. 

The next section describes the problem in more detail and states 
previous results. This is followed by a section presenting new 
results with a discussion of computational consideration and a section 
presenting an analytical example. 

2. Previous Results 


Object motion is usually modeled with a state space shaping filter 
where different classes of objects having different models may be con- 
sidered. It is assumed that the objects move independently of each 
other. Several sensors of different types may be in use with 
associated sensor error models. The parameters of the models in 
previous work have usually been assumed known. This paper includes 
unknown parameters entering the measurements linearly through 
initial conditions and forcing terms according to the following model 
for an object in a given class: 


>^k+l ■ \\ * ®2k‘>2 

^ -■ \ 

*0 ' * "or 

k = 0, 1 , . . . , M 


b = 



Fixed but unknown parameter vector 


(2-1) 
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W|^ - Zero mean, white, Gaussian process noise with covariance Qj^ 

V|^ - Zero mean, white, Gaussian measurement noise with co- 
variance R|^ 

Xor " Zero mean and Gaussian with covariance 

or 

For the remainder of this section - 0 , but b^ and bg are 

nonzero in the next section. 

The data association problem is displayed by Figure 2-1. Given 



-®— Object path with dot at measurement time t^ 
X Position measurement 
— Hypothesized data association 

Figure 2-1. Data Association Hypotheses 


172 



the data, there are a number of reasonable data associations. Many 
factors may have to be considered. Different sensors or the same 
sensor at different times may have a different coverage region, false 
alarm and detection errors may occur, and the number of objects may 
be initially unknown. In order to focus attention on new results, 
it is assumed in the remainder of the paper that the sensor coverage 
region includes all objects at each measurement time and that there 
are no false/alarm detection errors. Further, it is assumed that no 
a ■priori probabilistic information is available as to the correctness 
of a given data association. 

Previous results compute the statistical likelihood function with 
a bank of Kalman filters as shown in Figure (2-2) and maximize the 
likelihood over hypothesized data associations. The hypothesizer 
selects 



Figure 2-2. Likelihood Maximization 

the data association, and then the data is whitened by a bank of Kalman 
filters using appropriate object class models. Gonsequently, negative 
twice the log likelihood can be computed efficiently by the equation 
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-2£np{z;Hj) 




- I l 1 * C P.j 

•i,k ^i,k 


k 


'l,k > 


( 2 - 2 ) 


+ m tnZv 


” Residual and residual covariance for 
* hypothesis H. , object i , at time k 


m - Total number of scalar measurements. 


The data dependent part of equation (2-2) is chi-square with degrees 
of freedom m (x^) . Then the likelihood is maximized by minimizing 
equation (2-2). In this manner a statistically reasonable hypothesis 
is selected and tracking is performed by the bank of Kalman filters. 

3 . New Resul ts 


The contribution of this paper is the inclusion of fixed but un- 
known parameters in the object motion model. The generalized likeli- 
hood approach provides a solution by maximizing the likelihood over the 
parameters before maximizing over hypotheses. This is accomplished by 
replacing the Kalman filter for object i in Figure 2-2 by the equation 
shown in block diagram form in Figure 3-3 where the i and j sub- 
scripts are suppressed. The Kalman filter is based on the model of 
equation 2-1 with b = 0 This is indicated by subsripting filter computed 

quantities- with "fb - 0" The filter residuals and residual 
covariance are used in a set of auxilTiary equations to compute a maxi- 
mum likelihood estimate b, with estimate error covariance Pc 

k b|j 

These results ane used in turn to compute the portion of the likelihood 
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Figure 3-3, fiene^'aUzed Likelihood Computation 


for the object under consideration. Then using the filter estimate 
for the state and the estimate for b recovered from the filter 
residuals, the minimum variance and unbiased estimate Xj^ of the state 
and its estimate error covariance P- are computed. 

The estimate for b and the estimate uncertainty are given by 


< I "I I "I ’’i' 

fb=0 fb=0 


fb=0 


Tu ^ C0„^n, \h<^ 

^10 = i°nxn ] 

I c 

Vi 

- Kalman gain at time 
n^jOg - Number of elements in b^ and bg 
n - Number of elements in x. 


(3-1) 
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Note that the above equations can be implemented recursively. The 
minimized negative twice log likelihood for the given object is 
expressed as 


-2iin p(z;b|^) tn|P~ I + 

^ fb=0 fb=0 


■k 
fb=0 


) 


- bJi + m tn 2 tt 


(3-2) 


Now, including the i and j' subscripts, the generalized negative 
twice log likelihood over all objects is 


•2 in 

= Id I + z{k ^ik ^ 


^ ^ fb=0 fb=0 

fb=0 ^ fb=0 ^ 


i M h'i 


M 


iM 


(3-3) 


+ m tn 2 tt 


The data association is accomplished by minimizing the above equation 
over Hj . It can be shown that the data dependent part of equation 
3-3 is chi-square with degrees of freedom equal the total number of 
scalar measurements minus the total number of elements in all the b^ 
and bg vectors under hypothesis Hj . 

The generalized likelihood has been computed but valid state 
estimates for the objects still remain to be specified. It can be 
shown that the following equations give the minimum variance unbiased 
estimate and its uncertainty for a given object, again suppressing the 
i, j subscripts: 
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It is worthwhile to compare the new results of this section with 
the previous results of Section 2 where the noninformative prior 
approach is used to deal with the parameters b . The distribution of 
the data dependent part of the negative twice log likelihood is chi- 
square in the total number of scalar measurements for the previous 
results but is chi-square in the number of scalar measurements minus 
the total number of elements in all the b^ and b 2 vectors for the 
new results. However, it can be shown that the actual numerical value 
of the data dependent parts is the same. On the other hand, compari- 
son of equations (3-3) and (2-2) shows that the model dependent parts 
are different. Further, it can be shown that the state estimate 
obtained from the Kalman filter in the previous results is the same as 
that given by equation (3-4). There are clearly basic similarities 
and basic differences between the two approaches and caution must be 
exercised in selecting which approach is most appropriate. 

Numerical considerations are of major importance to the application 
of the new results as they are to previous results. It is easily 
verified that even for a modest number of objects and measurement times 
that a combinatoric explosion occurs in the number of hypothesized 
data associations. Much use can be made of pre-processing techniques 
to greatly limit the number of hypotheses. For example, if two posi- 
tion measurements require an object to exceed its maximum possible 
velocity then those measurements can not both come from the same object. - 
Also examination of filter residuals directly can be used to exclude 
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clearly unlikely data associations. In general, however, such pre- 
processing does not appear to stop the combinatoric explosion, but only 
to slow it. Thus, suboptimal approaches appear to be required for 
most realistic problems. For example, data association can be performed 
at each measurement time saving only tracks from the "best" hypotheses 
up to a maximum number of tracks. It is more efficient to save tracks 
than hypotheses since two hypotheses may have many tracks in conroon. 
Another technique is to again perform data association at each point in 
time but consider only hypotheses that can be generated by back track- 
ing N steps in time. Some mathematical progranming techniques are 
also applicable such as integer programming techniques (Reference 3). 
Ultimately, the best technique is determined by the details of the 
specific problem at hand. In general, what is desired as an algorithm 
with storage requirements that do not grow in time and with execution 
time that grows only linearly in time. 

4. Analytical Example 

A simple example is presented to provide some "feeling" for the 
generalized likelihood approach. The example is actually somewhat 
more general than the theory of the preceding section in that a power 
spectral density parameter is to be estimated from the data. 

The objects under observation move in a plane as random walks 
about straight mean paths with constant mean velocity where parameters 
describing the mean path are fixed but unknown. The power spectral 
density of the process noise will be considered known in one case and 
unknown in another. Measurements are taken at a constant interval 
At and are accurate enough to be considered noise free. It is 
assumed that all objects are observed at each measurement time and 
that no false alarm or detection errors occur. 

The analysis model for the i^ ■ object under the j ■ hypothesis 
is 
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k = 

r, . r„ . , Us u, . - Fixed but unknown parameters 

I 5 1 > 0 U > 1 9 0 V9 ^ 9 1 

w, . . , w„ > . - zero mean, white Gaussian each having variance 

I j J 5 K £ , I 9 K 

At and uncorrelated with each other 

The first measurement estimates initial position exactly so it will no 
longer be considered unknown. Using Kalman filter equations, the nega- 
tive twice log likelihood for the object of equation (4-1) is found to 
be 


-2 an p !,»••• 5 ., Q.) + constant 


« » / 
r ^^1 i k“^l i k-r^l i^^^' 

I 2 m Q^. At + ■ ' ' — LiJ 


M 

I 

k=l 


Q^.At 


(4-2) 


+ H,i,k'^2.i,k-T^2,i^^^‘ 

Q^At 


Now minimize equation (4-2) with respect to u^ ^ and ^ to obtain 
the maximum likelihood estimates 
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''J 1 In j j 

" MAt ■ ^l.i.k -1 

1 M . 

GJ = J_ y 2 ^ - 7J 

2 ,i MAt 2 ,i,k ^ 2 ,i,k-l 


(4-3) 


These estimates are valid whether Q^. is known or not. Now assume 
is unknown and minimize equation (4-2) to obtain 


M 


^i "2MAt ^^l.i.k'^l.i.k-r'^l.i^^^ 




(4-4) 


Now consider the case where the Q- are known and form the general- 
ized negative twice log likelihood over all objects minimized over u^ ^ 
and u« . 

4 j 1 


-2 tn p (z; 


J_ 

At 




•*,Ui .,u^ ^»***) + constant 

I j 1 4 9 I 


M 


^^1 ,i ,k“^l ,i ,k-f M ^^1 ,i ,r^l ,i ,1-1 ^ 


^4,i ,k'4,i ,k-f M ^4,i ,r4,i ,i-l 


(4-5) 


For the case of Q. unknown, the generalized likelihood is maximized 
by minimizing the following expression 


180 



(4-6) 


^ ^2 , i . k"^2 , i , k-f FlJ ^ ^ i , a’^ 2 , i , a-I ^ ^ ^ 


Equations (4-5) and (4-6) have a term that sums over time in common. 
This common term is simply a measure of the variation of position 
change about mean position change for a given object. Both equations 
qualitatively attempt to perform data association by minimizing the 
variation in position change over all objects. However, In, equation 
(4-5) where Q> is known, the variation in position change is weighted 
heavily for small Q. . This is intui£ively reasonable since a small 
variation would be expected. It might seem reasonable that if Q,. 
is unknown that equation (4-5) should still apply but with 1 . 
However, the variation in position change for objects that truly have a 
small would be swamped out by that for objects with large Q^. . 

The multiplication of equation (4-6) keeps this swamping effect from 
occurring. 
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